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Abstract. How and when did the first generation of stars form at the end of the cosmic dark 
ages? Quite generically, within variants of the cold dark matter model of cosmological structure 
formation, the first sources of light are expected to form in ~ 1O 6 M0 dark matter potential 
wells at redshifts z ^ 20. I discuss the physical processes that govern the formation of the first 
stars. These so-called Population III stars are predicted to be predominantly very massive, and 
to have contributed significantly to the early reionization of the intergalactic medium. Such 
an early reionization epoch is inferred from the recent measurement of the Thomson optical 
depth by the WMAP satellite. I address the importance of heavy elements in bringing about 
the transition from an early star formation mode dominated by massive stars, to the familiar 
mode dominated by low mass stars, at later times, and present possible observational probes. 
This transition could have been gradual, giving rise to an intermediate-mass population of still 
virtually metal-free stars ("Population II. 5"). These stars could have given rise to the peculiar 
class of black-hole forming supernovae inferred from the abundance pattern of extremely iron- 
poor stars. 
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1. Introduction 

One of the grand challenges in modern cosmology is posed by the question: How did the 
first stars in the universe form, what were their physical properties, and what was their 
impact on cosmic history (e.g., Bromm & Larson 2004)? The first stars, formed at the end 
of the cosmic dark ages, ionized (e.g., Wyithe & Loeb 2003; Cen 2003) and metal-enriched 
(e.g., Furlanetto & Loeb 2003) the intergalactic medium (IGM) and consequently had 
important effects on subsequent galaxy formation (e.g., Barkana & Loeb 2001) and on the 
large-scale polarization anisotropics of the cosmic microwave background (Kaplinghat et 
al. 2002). When did the cosmic dark ages end? In the context of popular cold dark matter 
(CDM) models of hierarchical structure formation, the first stars are predicted to have 
formed in dark matter halos of mass ~ 10 6 Af Q that collapsed at redshifts z ~ 20 — 30 
(e.g., Barkana & Loeb 2001; Yoshida et al. 2003). 

Results from recent numerical simulations of the collapse and fragmentation of pri- 
mordial clouds suggest that the first stars were predominantly very massive, with typical 
masses M* ^ 1OOM (Bromm, Coppi, & Larson 1999, 2002; Nakamura & Umemura 
2001; Abel, Bryan, & Norman 2002). Despite the progress already made, many impor- 
tant questions remain unanswered. An example for an open question is: How does the 
primordial initial mass function (IMF) look like? Having constrained the characteristic 
mass scale, still leaves undetermined the overall range of stellar masses and the power- 
law slope which is likely to be a function of mass. In addition, it is presently unknown 
whether binaries or, more generally, clusters of zero-metallicity stars, can form. What 
is the nature of the feedback that the first stars exert on their surroundings? The first 
stars are expected to produce copious amounts of UV photons and to possibly explode 
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Figure 1. Collapse and fragmentation of a primordial cloud. Shown is the projected gas density 
at a redshift z ~ 21.5, briefly after gravitational runaway collapse has commenced in the center 
of the cloud. Left: The coarse-grained morphology in a box with linear physical size of 23.5 pc. 
Right: The fine-grain morphology in a box with linear physical size of 0.5 pc. The central density 
peak, vigorously gaining mass by accretion, is accompanied by a secondary clump. (Adapted 
from Bromm & Loeb 2004.) 

as energetic hypernovae. These negative feedback effects could suppress star formation 
in neighboring high-density clumps. 

Predicting the properties of the first sources of light, in particular their expected lu- 
minosities and spectral energy distributions, is important for the design of upcoming 
instruments, such as the James Webb Space Telescope (JWST) f, or the next generation 
of large (> 10m) ground-based telescopes. The hope is that over the upcoming decade, 
it will become possible to confront current theoretical predictions about the properties 
of the first stars with direct observational data. 

2. Population III star formation 

The metal-rich chemistry, magnetohydrodynamics, and radiative transfer involved in 
present-day star formation is complex, and we still lack a comprehensive theoretical 
framework that predicts the IMF from first principles (see Larson 2003 for a recent 
review). Star formation in the high redshift universe, on the other hand, poses a theo- 
retically more tractable problem due to a number of simplifying features, such as: (i) the 
initial absence of heavy metals and therefore of dust; and (ii) the absence of dynamically- 
significant magnetic fields, in the pristine gas left over from the big bang. The cooling 
of the primordial gas does then only depend on hydrogen in its atomic and molecu- 
lar form. Whereas the initial state of the star forming cloud is poorly constrained in 
the present-day interstellar medium, the corresponding initial conditions for primordial 
star formation are simple, given by the popular ACDM model of cosmological structure 
formation. 

How did the first stars form? A complete answer to this question would entail a the- 
oretical prediction for the Population III IMF, which is rather challenging. Let us start 
by addressing the simpler problem of estimating the characteristic mass scale of the first 
stars. This mass scale is observed to be ~ 1M Q in the present-day universe. To investigate 

f See http:// ngst.gsfc.nasa.gov. 
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the collapse and fragmentation of primordial gas, we have carried out numerical simula- 
tions, using the smoothed particle hydrodynamics (SPH) method. We have included the 
chemistry and cooling physics relevant for the evolution of metal-free gas (see Bromm et 
al. 2002 for details). Improving on earlier work (Bromm et al. 1999, 2002) by initializing 
our simulation according to the ACDM model, we focus here on an isolated overdense 
region that corresponds to a 3a— peak: a halo containing a total mass of 1O 6 M0, and 
collapsing at a redshift z v ; r ~ 20 (Bromm & Loeb 2004). In Figure 1 (left panel), we show 
the gas density within the central ~ 25 pc, briefly after the first high-density clump has 
formed result of gravitational runaway collapse. 

How massive were the first stars? Star formation typically proceeds from the 'inside- 
out', through the accretion of gas onto a central hydrostatic core. Whereas the initial 
mass of the hydrostatic core is very similar for primordial and present-day star formation 
(Omukai & Nishi 1998), the accretion process - ultimately responsible for setting the final 
stellar mass, is expected to be rather different. On dimensional grounds, the accretion 
rate is simply related to the sound speed cubed over Newton's constant (or equivalently 
given by the ratio of the Jeans mass and the free-fall time): M acc ~ cf/G a T 3 / 2 . A 
simple comparison of the temperatures in present-day star forming regions (T ~ 10 K) 
with those in primordial ones (T ~ 200 — 300 K) already indicates a difference in the 
accretion rate of more than two orders of magnitude. 

Our high-resolution simulation enables us to study the three-dimensional accretion 
flow around the protostar (see also Omukai & Palla 2001, 2003; Ripamonti et al. 2002; 
Tan & McKee 2003). We allow the gas to reach densities of 10 12 cm~ 3 before being 
incorporated into a central sink particle. At these high densities, three-body reactions 
(Palla, Salpcter, & Stahler 1983) have converted the gas into a fully molecular form. In 
Figure 2, we show how the molecular core grows in mass over the first ~ 10 4 yr after its 
formation. The accretion rate (left panel) is initially very high, M acc ~ 0.1M Q yr -1 , and 
subsequently declines with time. The mass of the molecular core (right panel), taken as 





time [yr] 

Figure 2. Accretion onto a primordial protostar. The morphology of this accretion flow is shown 
in Fig. 1. Left: Accretion rate (in Mq yr -1 ) vs. time (in yr) since molecular core formation. 
Right: Mass of the central core (in Mq) vs. time. Solid line: Accretion history approximated as: 
M* oc t ' 75 at t < 10 3 yr, and M» oc t 0A afterwards. Using this analytical approximation, we 
extrapolate that the protostellar mass has grown to ~ 12OM after ~ 10 5 yr, and to ~ 500M Q 
after ~ 3 x 10 6 yr, the total lifetime of a very massive star. (Adapted from Bromm & Loeb 
2004.) 



4 



V. Bromm 



an estimator of the proto-stellar mass, grows approximately as: M* ~ J M acc dt cx t 75 at 
t < 10 3 yr, and M* cx t 0A afterwards. A rough upper limit for the final mass of the star is 
then: M,(t=3x 10 6 yr) ~ 5OOM . In deriving this upper bound, we have conservatively 
assumed that accretion cannot go on for longer than the total lifetime of a very massive 
star (VMS). 

Can a Population III star ever reach this asymptotic mass limit? The answer to this 
question is not yet known with any certainty, and it depends on whether the accretion 
from a dust-free envelope is eventually terminated by feedback from the star (e.g., Omukai 
& Palla 2001, 2003; Omukai & Inutsuka 2002; Ripamonti et al. 2002; Tan & McKee 
2003). The standard mechanism by which accretion may be terminated in metal-rich gas, 
namely radiation pressure on dust grains (e.g., Wolfire & Cassinelli 1987), is evidently 
not effective for gas with a primordial composition. Recently, it has been speculated that 
accretion could instead be turned off through the formation of an H II region (Omukai & 
Inutsuka 2002), or through the radiation pressure exerted by trapped Lya photons (Tan 
& McKee 2003). The termination of the accretion process defines the current unsolved 
frontier in studies of Population III star formation. 

3. Second generation stars 

How and when did the transition take place from the early formation of massive stars 
to that of low-mass stars at later times? In contrast to the formation mode of massive 
stars (Population III) at high rcdshifts, fragmentation is observed to favor stars below a 
solar mass (Population I and II) in the present-day universe. The transition between these 
fundamental modes is expected to be mainly driven by the progressive enrichment of the 
cosmic gas with heavy elements, which enables the gas to cool to lower temperatures. 
The concept of a 'critical mctallicity', Z cr j t , has been used to characterize the transition 
between Population III and Population II formation modes, where Z denotes the mass 
fraction contributed by all heavy elements (Omukai 2000; Bromm ct al. 2001; Schneider 
et al. 2002: Schneider et al. 2003; Mackey, Bromm, & Hernquist 2003; Yoshida, Bromm, 
& Hernquist 2004). These studies have constrained this important parameter to only 
within a few orders of magnitude, Z crlt ~ 10~ 6 — 10~ 3 Zq, under the implicit assumption 
of solar relative abundances of metals. This assumption is likely to be violated by the 
metal yields of the first SNe at high-rcdshifts, for which strong deviations from solar 
abundance ratios are predicted (e.g., Heger & Woosley 2002; Qian & Wasserburg 2002; 
Umeda & Nomoto 2002, 2003). 

Recently, we have shown that the transition between the above star formation modes 
is driven primarily by fine-structure line cooling of singly-ionized carbon or neutral 
atomic oxygen (Bromm & Loeb 2003). Earlier estimates of Z cr i t which did not explic- 
itly distinguish between different coolants arc refined by introducing separate critical 
abundances for carbon and oxygen, [C/H] cr it and [0/H] cr i t , respectively, where [A/H] = 
\og w (NA/Nu) — log 10 (A r A/A^H)o- Since C and O are also the most important coolants 
throughout most of the cool atomic ISM in present-day galaxies, it is not implausible 
that these species might be responsible for the global shift in the star formation mode. 
Under the temperature and density conditions that characterize Population III star for- 
mation, the fine-structure lines of O I and C II dominate over all other metal transitions. 
Cooling due to molecules becomes important only at lower temperatures, and cooling 
due to dust grains only at higher densities (e.g., Omukai 2000, Schneider et al. 2003). 
Numerically, the critical C and O abundances are estimated to be: [C/H] cr it — — 3.5±0.1 
and [0/H] crit ~ -3.1 ±0.2. 

Even if sufficient C or O atoms are present to further cool the gas, there will be a 
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Figure 3. Evolution of HD abundance in four distinct cases. Solid line: Gas compressed and 
heated by a SN explosion with u B h~ 100 km s _1 . Dotted line: Gas shocked in the build-up of a 
3<r fluctuation dark matter halo collapsing at z ~ 15. Dash-dotted line: Gas collapsing inside a 
relic HII region, which is left behind after the death of a very massive Pop III star. Dashed line: 
Gas collapsing inside a minihalo at z ~ 20. In this case, the gas does not experience a strong 
shock, and is never ionized. Contrary to the other three cases, where the gas went through a fully 
ionized phase, HD cooling is not important here. The critical HD abundance, shown by the bold 
dashed line, is defined such that primordial gas is able to cool to the CMB temperature within 
the fraction of a Hubble time. The CMB sets a minimum floor to the gas temperature, because 
radiative cooling below this fllor is thermodynamically not possible. The HD abundance exceeds 
the critical value in a time which is short compared to the Hubble time for all fully-ionized, 
strongly shocked cases. (Adapted from Johnson & Bromm 2005.) 

minimum attainable temperature that is set by the interaction of the atoms with the 
thermal CMB: Tomb = 2.7 K(l + z) (e.g., , ). At z ~ 15, this results in a characteristic 
stellar mass of M* - 2OM (n / /lO 4 cm- 3 )- 1 /2 j wliere nf > 1Q 4 

cm 3 is the density at 

which opacity prevents further fragmentation. It is possible that the transition from the 
high-mass to the low-mass star formation mode was modulated by the CMB tempera- 
ture and was therefore gradual, involving intermediate-mass ('Population II. 5') stars at 
intermediate redshifts (). This transitional population could give rise to the faint SNe 
that have been proposed to explain the observed abundance patterns in metal-poor stars 
(Umeda & Nomoto 2002, 2003). When and how uniformly the transition in the cosmic 
star formation mode did take place was governed by the detailed enrichment history of 
the IGM. This in turn was determined by the hydrodynamical transport and mixing of 
metals from the first SN explosions (e.g., Mori, Ferrara, & Madau 2002; Bromm, Yoshida, 
& Hernquist 2003; Scannapieco, Schneider, & Ferrara 2003; Wada & Venkatesan 2003). 

Recently, the additional boost to the cooling of still metal-free gas provided by HD has 
been investigated (e.g., Johnson & Bromm 2005, and references therein). If the primordial 
gas goes through a strongly shocked, fully ionized phase prior to the onset of protostellar 
collapse, cooling is possible down to the temperature of the CMB which sets the minimum 
floor accessible via radiative cooling (see Figure 3) . The lower temperatures in turn could 
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allow the fragmentation into intermediate-mass stars, with masses of order a few tens of 
Mq, giving rise to a possible "Population II. 5" (see Figure 4). 



4. Stellar fossils 

It has long been realized that the most metal-poor stars found in our cosmic neighbor- 
hood would encode the signature from the first stars within their elemental abundance 
pattern. For many decades, however, the observational search has failed to discover a 
truly first-generation star with zero metallicity. Indeed, there seemed to have been an 
observational lower limit of [Fc/H] ~ —4 (e.g., Carr 1987). In view of the recent theoreti- 
cal prediction that most Population III stars were very massive, with associated lifetimes 
of ~ 10 6 yr, the failure to find any 'living' Population III star in the Galaxy is not sur- 
prising, as they would all have died a long time ago. Furthermore, theory has predicted 
that star formation out of extremely low-metallicity gas, with Z < Z cr ;t ~ 10~ 35 Z Q , 
would be essentially equivalent to that out of truly primordial gas. Again, this theoreti- 
cal prediction was in accordance with the apparent observed lower-metallicity cutoff. 

Recently, however, this simple picture has been challenged by the discovery of the star 
HE0107-5240 with a mass of 0.8M© and an iron abundance of [Fe/H] = -5.3 (Christlieb 
et al. 2002; Frebel et al. 2005). This finding indicates that at least some low mass stars 
could have formed out of extremely low-metallicity gas. Does the existence of this star 
invalidate the theory of a metallicity threshold for enabling low-mass star formation? A 
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Figure 4. Characteristic mass of stars as a function of redshift. Pop III stars, formed out of 
metal-free gas, and not going through a fully ionized phase prior to the onset of collapse, have 
typical masses of M» ~ IOOMq. Pop II stars, formed out of already metal-enriched gas, are 
formed at lower redshifts, and have typical masses of M* ~ 1M . Pop II. 5 stars, formed out 
of strongly shocked, but still virtually metal-free gas, are hypothesized to have typical masses 
that are intermediate between Pop III and Pop II with M* ~ IOMq. (Adapted from Johnson & 
Bromm 2005.) 
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Figure 5. Observed abundances in low-metallicity Galactic halo stars. For both carbon (upper 
panel) and oxygen (lower panel), filled circles correspond to samples of dwarf and subgiant 
stars, and open squares to a sample of giant stars (see for references) . The dashed lines indicate 
the predicted critical carbon and oxygen abundances. Highlighted (cross) is the location of the 
extremely iron-poor giant star HE0107-5240. (Adapted from Bromm & Loeb 2003.) 



possible explanation (Umeda & Nomoto 2003) could lie in the unusually high abundances 
of carbon and oxygen in HE0107-5240. 

In Figure 5, the theoretical C and O thresholds (Bromm & Loeb 2003) are compared 
to the observed abundances in metal-poor dwarf and giant stars in the halo of our Galaxy 
(see Bromm & Loeb 2003 for references). As can be seen, all data points lie above the 
critical O abundance but a few cases lie below the critical C threshold. All of these low 
mass stars are consistent with the model since the corresponding O abundances lie above 
the predicted threshold. The sub-critical [C/H] abundances could have either originated 
in the progenitor cloud or from the mixing of CNO-processed material (with carbon 
converted into nitrogen) into the stellar atmosphere during the red giant phase. Note 
that the extremely iron-poor star HE0107-5240 has C and O abundances that both lie 
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above the respective critical levels. The formation of this low mass star (~ 0.8M Q ) is 
therefore consistent with the theoretical framework considered by Bromm & Loeb (2003). 

The lessons from stellar archaeology on the nature of the first stars are likely to increase 
in importance, since greatly improved, large surveys of metal-poor Galactic halo stars 
are under way, or are currently being planned. 
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